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Abstract 

We study the IIB matrix model, which is conjectured to be a nonperturbative defini- 
tion of superstring theory, by introducing an integer deformation parameter u which 
couples to the imaginary part of the effective action induced by fermions. The de- 
formed IIB matrix model continues to be well-defined for arbitrary u, and it preserves 
gauge invariance, Lorentz invariance, and the cluster property. We study the model 
at v = oo using a saddle-point analysis, and show that ten-dimensional Lorentz in- 
variance is spontaneously broken at least down to an eight-dimensional one. We argue 
that it is likely that the remaining eight-dimensional Lorentz invariance is further bro- 
ken, which can be checked by integrating over the saddle-point configurations using 
standard Monte Carlo simulation. 
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Recent developments in string dualities have suggested that the known five types of su- 
perstring theories in ten dimensions and M-theory in eleven dimensions are just different 
microscopic descriptions of the same underlying physics. This means that if any of the 
descriptions can be defined nonperturbatively, it is as good as any other to study the dy- 
namics of the universality class of string/M theories. A celebrated first step towards this 
end has been made by Ref. |L], where the so-called Matrix theory has been proposed as 
a nonperturbative definition of M-theory in the infinite momentum frame. Soon after, the 
IIB matrix model @ has been proposed as a nonperturbative definition of type IIB super- 
string theory. Although there are many evidences that support these conjectures, there is 
no direct proof that these models reproduce the string/M theories perturbativelyQ. If these 
conjectures are true, there is a hope to understand all the fundamental questions about the 
Standard Model, including the space-time dimensionality, gauge group, matter contents, the 
hierarchy problem, the cosmological constant problem, and so on, in terms of the dynamics 
of these models. 

So far, the understanding of the vacuum of these models is quite limited. In this Letter, 
we attempt to extract information about the dimensionality of the space-time which is 
dynamically generated in the IIB matrix model. The IIB matrix model is a supersymmetric 
matrix model obtained formally by taking a zero- volume limit [|J of ten-dimensional SU(iV) 
super Yang-Mills theory. Unlike in field theories, the space-time is treated in this model 
as a dynamical object, which is represented by ten bosonic N x N hermitian matrices, 
where N should be sent to infinity eventually. Therefore, the model has a potential to 
explain the dynamical origin of the dimensionality of our space-time, which, in the Standard 
Model, has to be given as an input parameter. If we are ever going to explain our four- 
dimensional space-time in this way, the vacuum of the IIB matrix model should be dominated 
by configurations which have only four-dimensional extent^. This means, in particular, that 
the ten-dimensional Lorentz invariance of the model should be spontaneously broken to a 
four-dimensional one.f] 

Such an issue can in principle be addressed by performing Monte Carlo simulation, since 
the model is completely well defined for arbitrary iV without any cutoff |7], |8|, |9|. Indeed 
the large-iV dynamics of the four- dimensional version of the IIB matrix model, which can 
be obtained by a zero-volume limit of four- dimensional super Yang-Mills theory, has been 
understood to a considerable extent through Monte Carlo simulation |TI]] . Comparison of the 
1 See, however, Ref. || as an attempt in this direction. 

2 This is similar in spirit to a description of our four-dimensional space-time as a "brane" in a higher- 
dimensional non-compact space-time, which is proposed Q as an alternative to a more conventional com- 
pactification approach. 

3 The possibility of a non-trivial vacuum structure in string theory has been considered also in [0. 
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results with those for the corresponding bosonic theory [TT[, revealed that supersymmetry 
indeed plays an important role in the dynamics of large- A reduced models as a nonpertur- 
bative definition of string theories. Unfortunately, such a direct approach does not work 
in the IIB matrix model as it stands, due to the fact that the effective action induced by 
fermions is generically complex, whereas in the above mentioned four-dimensional version, 
it is real. In general, when the action of a theory has a non-zero imaginary part, the number 
of configurations needed to extract any information increases as exponential of the system 
size. This is the notorious sign problem, which occurs also in many other interesting systems 
related to particle physics, such as theories with a chiral fermion (as in the present case), 
theories with a 9 vacuum, Chern-Simons theories and theories with a finite baryon number 
density. We emphasize, however, that the problem is a purely technical one, and indeed its 
complete solution in some particular class of systems has been obtained recently f[2 |. 

In order to search for a way out of this difficulty, let us consider a deformation of the 
IIB matrix model by introducing an integer parameter z/, which couples to the imaginary 
part of the effective action induced by fermions. The deformed model, which reduces to 
the IIB matrix model at v = 1, is well-defined for arbitrary v. It preserves both gauge 
invariance and Lorentz invariance, and moreover, it preserves the cluster property, which 
is an important consequence of the supersymmetry of the original model. At v — 0, the 
sign problem disappears and standard Monte Carlo techniques become applicable. Monte 
Carlo studies up to N = 512 show, however, that the ten-dimensional Lorentz invariance 



is not spontaneously broken [13| , which suggests that if Lorentz invariance is spontaneously 
broken in the original model, the imaginary part of the effective action must play a crucial 
role. This motivated us to consider the opposite extreme, namely v = oo. 

At v = oo, the integration over the bosonic matrices is dominated by the saddle-points 
of the imaginary part of the effective action. We find that the saddle-points are given 
by configurations which have only eight-dimensional extent. This implies that the ten- 
dimensional Lorentz invariance is spontaneously broken at least to an eight-dimensional 
one. Since we still have to integrate over the saddle-point configurations, the question arises 
whether the remaining eight-dimensional Lorentz invariance further breaks down, say, to a 
four- dimensional one. To this end, we study the hessian at configurations which have only 
(i-dimensional extent (d < 8) and find that it is zero of order (8 — d) {2 (N 2 — 1) — 16}. In 
other words, the imaginary part of the effective action becomes more stationary for smaller 
d. This gives a huge enhancement to configurations with smaller d, which is found to cancel 
exactly the entropical barrier against having such configurations. Considering the effect of 
the real part of the effective action, which is studied in Ref. f| using perturbation theory, 
it is plausible that the vacuum of the model at v = oo is actually given by d < 8. We also 
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find a lower bound on d as d > 2. Further information about the vacuum can be obtained 
by performing explicitly the integration over the saddle-point configurations by Monte Carlo 
simulation. Naively, one might consider that the sign problem, which already exists in the 
original model {y — 1), becomes maximally severe at v — oo. This is not the case, as we will 
see. 

The IIB matrix model is formally a zero-volume limit of ten-dimensional pure M = 1 
supersymmetric Yang-Mills theory. The action, therefore, is given by 

S = St + Sf , 
S b = -\tx {[A^A v f) , 

S f = ~tx {^{CYX^M) ■ (i) 

(fj, = 1, ■ • • , 10) and ip a (a = 1, ■ ■ ■ , 16) are N x N traceless hermitian matrices, which 
can be expanded in terms of the generators t a of SU(iV) as 

N 2 -l N 2 ~l 
a=l a=l 

where A ^ is a real variable and ip^ is a real Grassmann variable. We have made a Wick 
rotation in the action (|l|), so that the metric has Euclidean signature. The 16 x 16 matrices 
are ten-dimensional gamma matrices after Weyl projection, and the unitary matrix C is 
a charge conjugation matrix satisfying 

CT^ = (T,) T ; C T = C. (3) 

The model has a manifest ten-dimensional Lorentz invariance, by which we actually mean 
an SO(10) invariance. A^ transforms as a vector and ip a transforms as a Majorana-Weyl 
spinor. 

Pure M = 1 supersymmetric Yang-Mills theory can be also defined in 3D, 4D and 6D, 
as well as in 10D. Hence, by taking a zero-volume limit of these theories, we can define 
supersymmetric large- N reduced models, which are D = 3,4,6 versions of the IIB matrix 
model. A nontrivial question then concerns whether the integration over the bosonic ma- 
trices is convergent, since the integration domain for hermitian matrices is non-compact. A 
potential danger of divergence, even for finite N, exists when the eigenvalues of A^ become 
large. This issue has been addressed in Ref. || using a one-loop perturbative argument. 
When all the eigenvalues are well separated from each other, one can expand the matrices 
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and ip a around diagonal matrices as 

where a M and y> Q are matrices containing only off-diagonal elements. One can then integrate 
over the off-diagonal elements a M and <p a up to one-loop, by fixing the gauge properly and 
including the Faddeev-Popov ghosts. The integration over the bosonic off-diagonal elements 
dp and the Faddeev-Popov ghosts, gives a logarithmic attractive potential between all the 



pairs of x% This potential, however, is exactly cancelled by the contribution of the 

fermionic off-diagonal elements (p a due to supersymmetry. This cancellation is responsible 
for the cluster property of the model 0], which is important for the interpretation of the 
model as a string theory. In order to calculate the effective potential for x^, one still has to 
integrate over the fermionic diagonal elements £j a , which is nontrivial. In Ref. ||, it has been 
shown that the effective potential for x^ can be given by a branched-polymer like attractive 
interaction among N points in D-dimensional space-time represented by x^ (i — 1, • • • , N). 
Now counting the power of x« in the partition function, the integration measure for x^ gives 
D (N — 1) and the integration over £ ia gives — %p (N — 1), where p = 2 (D — 2) is the number 
of real components of the spinor considered for each D. Therefore, naively one concludes 
that when 

D (N - 1) - ~p(N - 1) = 2 (3 - D) (N - 1) > , (5) 

the integration over Xi^ is divergent and otherwise convergent. This means that the integral 
is divergent for D = 3 and convergent for D = 4,6, 10, irrespectively of N. This conclusion 
is in agreement with an exact result available for # = 2 @ and a numerical result obtained 
for JV = 3 H. For D = 4, Monte Carlo simulations show that the conclusion extends to the 



large- N limit |T0|. Therefore, it is conceivable that the above conclusion obtained by the 
one-loop argument holds in general, although it is not a rigorous proof. 

Going back to the definition of the model fll]), let us first integrate over the fermionic 
matrices ip a , which induces an effective action for A^. When fermions are real Grassmann 
variables as in the present D = 10 case, the fermion integral yields a pfaffian of a p (N 2 — 1) x 
p (N 2 — 1) matrix. On the other hand, when fermions are complex Grassmann variables as in 
the D = 4, 6 case, the fermion integral yields a determinant of a p (N 2 — l)/2 x p (N 2 — 1) /2 
matrix. For D = 4, the fermion determinant is real positive, which allows a direct Monte 
Carlo study of the model [ID|. For D = 6, 10, the fermion determinant (D = 6) or pfaffian 



(D = 10) is complex in general, since fermions are essentially chiral in these cases. This 

causes the notorious sign problem. There are some exceptions^ when iV is small 0, [| . For 

4 This is because, when (iV 2 — 1) is smaller than D, one can always make a D-dimensional rotation to set 
= for fi = 1, • • • , D — (N 2 — 1), since one can regard as (N 2 — 1), D-dimensional real vectors. The 
statements in D = 10 for N — 2 and N — 3, for example, follow from ( |l4| ) and (|l^), respectively. 
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D — 6, the fermion determinant is real positive for A = 2. For D = 10, the pfaffian is real 
positive for iV = 2 and real (but not positive definite) for A = 3. Although we consider only 
the D — 10 case in what follows, the analysis can be readily applied to the D — 6 case as 
well. We also restrict ourselves to iV > 4 so that the pfaffian is complex generically. 
In D = 10, the fermion integral yields 

Z f [A] = J d^j e~ s f =PfM, (6) 

where 

M aa , bp = -i f ahc ( C r M ) Q/3 A£ (7) 

is a 16 (A 2 — 1) x 16 (iV 2 — 1) anti-symmetric matrix, regarding each of (aa) and (bj3) as 
a single index. The real totally-antisymmetric tensor f a b c gives the structure constants of 
SU(iV). In what follows, it proves convenient to work with an explicit representation of the 
gamma matrices given by 

T l = i cr 2 ® o- 2 <g> <7 2 ® cr 2 ; r 2 = i 2 <g> 2 ® 1 ® 0i ; r 3 = i 2 ® 2 ® 1 <g> 3 ; 

r 4 = % 2 ® 0i (g) 2 ® 1 ; r 5 = % 2 ® 3 ® 2 ® 1 ; r 6 = % a 2 ® 1 ® (J\ ® a 2 ; 

r 7 = i0 2 (g)l(g)03(g)0 2 ; r 8 = 2 0i(g)l(g)l(g)l; r 9 = i0 3 (g)l(g)l(g)l; 

rio = 1 ® 1 ® 1 <8> 1 , (8) 

for which the charge conjugation matrix C becomes a unit matrix. 

Some important properties of the pfaffian Pf A4(A) are in order. First of all, it is invariant 
under a ten-dimensional Lorentz transformation, A'^ = A M „A„, where A is an SO(10) matrix. 
It is also invariant under a gauge transformation A'^ = gA u g J< , where g is an SU(A r ) matrix. 
Under a parity transformation, 

Af = -A 10 ; Af = Ai (for i = l, •••,9), (9) 

the pfaffian becomes complex conjugate, 

PiM(A p ) = {PtM(A)}* , (10) 

since M(A P ) = M(A)* . As a consequence, the pfaffian Pf M.{A) is real, when A w = 0. 
This also means that det M = (Pf A^) 2 is real positive, when A 10 = 0. 

We next consider the case with A 9 = A 10 = 0. In this case, the pfaffian is actually equal 
to the fermion determinant of eight-dimensional Weyl fermion which we denote as det A4®. 
M (8) is an 8 (A 2 - 1) x 8 (A 2 - 1) matrix defined as = -i f abc (F ^A^, where the 

eight-dimensional gamma-matrices after Weyl projection are given, for example, by 

ri = % 2 ® 2 ® 2 ; r 2 = % a 2 ® 1 <g> G\ ; r 3 = % 2 ® 1 ® 3 ; r 4 = % o\ ® a 2 <8> i ; 

f 5 = % 03 ® 02 ® 1 ; f 6 = Z 1 ® 0i ® 2 ; f 7 = % 1 ® 3 ® 2 J f 8 = 1 ® 1 ® 1 . (11) 
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With this representation, one finds that Ai^ is a real matrix, which means that detTVf'- 8 ^ is 
real. However, it is not positive definite, as we have checked numerically. In fact, this is the 
case even for A$ = 0. However, if we take Aj = As = 0, it becomes positive definite, since 
in this case, we have detAi^ = | det A / f ( - 6 * ) | 2 , where det Ai^ is the fermion determinant of 
six-dimensional Weyl fermion ||. 

Finally, when A 3 = A A = ■ ■ ■ = A 10 = 0, one finds that Pf M(A) =0 |J. This can be 
proved in the following way. First we note that (Pf M) 2 = det M = | det U\ 16 , where U is an 
(N 2 - 1) x (N 2 - 1) matrix defined as U ab = f abc X c , where X c = A\+iA c 2 . Since U ab X b = 0, 
the matrix U has a zero-eigenvalue, and therefore det U = 0. This completes the proof. 

Using Lorentz invariance of the pfaflian Pf At(A), we summarize some of its important 
properties as follows. Let us first define sets of "degenerate" configurations Qd as 

Cl d = { {A^} ; nfAp = for 3nf(i = 1, • • • , 10 - d) linearly independent} . (12) 

Obviously, Qi C f2 2 C ■ • • C ^9. Then the statements are 

(a) When {A^} G fig , PiM(A) G K . (13) 

(b) When {A^} G fi 6 , PfM(A) > . (14) 

(c) When {A^} G fi 2 , Pf M(A) = . (15) 

Here we recall that in the IIB matrix model, the space-time is treated as a dynamical object 
represented by the bosonic matrices A^. In this regard, generic configurations in describe 
(^-dimensional space-time. It is very suggestive that the phase of the pfaflian is sensitive to 
the dimensionality d of the space-time. Indeed, the above results will play a crucial role in 
the following analysis. 



Let us write the effective action induced by fermion integral as r e ff = — ln(Pf Ai), which 
we decompose into the real part and the imaginary part rW as T c g = +iV^K Then 
we generalize the theory by introducing an integer parameter u, which couples only to 
the imaginary part of the effective action, as 

Z v = J dA e -s b [A]-r^[A] e -i„r«[A] (16) 

The parameter v has to be an integer since is defined only up to modulo 2ir. For v — 1, 
the model reduces to the original model. Since flips its sign under a parity transformation 
due to ([IUD, the models with v and —v are nothing but parity partners. Since and are 
both Lorentz invariant and gauge invariant separately, the generalized model is also both 
Lorentz invariant and gauge invariant. According to the one-loop perturbative argument 
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given below (Q), the deformed IIB matrix model continues to be well-defined for arbitrary 
v, and moreover, it preserves the cluster property, which is an important consequence of the 
supersymmetry of the original model. 

Another important point about the generalization we consider is that (Sb) can be calcu- 
lated exactly and the result is independent of the parameter v. For this, we recall that in 
the original supersymmetric large- N reduced models [D = 4, 6, 10), one obtains 

(S b ) = 1 -(d + ^j(N 2 -1)= 1 -(D- 1) (N 2 - 1) , (17) 

using a scaling argument similar to [fLl]| . Due to the fact that r^[A] is invariant under a 
scale transformation of A^, one can easily find that the above result is unaltered by the 
generalization to arbitrary v for D = 6, 10. In particular, this holds true even in the v — > oo 
limit. Note, however, that this does not necessarily imply that all the vacuum expectation 
values of the generalized model are independent of v. Even their convergence in the v — > oo 
limit is nontrivial. This point shall be clarified later. 

We prove a property of r^, which turns out to be essential in the analysis of the v — > oo 
limit. We first recall that the pfaffian Pf A4(A) is a polynomial of A® of order 8 (N 2 — 1). 
Hence, T^' is infinitely differentiable at a configuration, for which the matrix Ai is invertible. 
When the configuration belongs to Qd (d = 1, • • • , 8), we find that 

Qn p(i) l Qn 

r, Aai r, A a 2 ttt^ = ~o 7^A-a~ lmliadetM = for n = 1, • • • , (9 - d) . 

dAp\ dA£ ■ ■ ■ dAli 2 dA^\ dA£ ■ ■ ■ dA a ^ n 

(18) 

This is because, up to (9 — d)-th order of perturbations, the configuration stays within Qg, 
and therefore det Ai remains to be real positive due to fll3|) . 

In the v — > oo limit, the integration over A^ is dominated by the configurations which 
satisfy the saddle-point equation given by 

= ^ = -^ Imlndet - M = -^ ImTr {M- l B^) , (19) 
where B^ a) is an A-independent 16 (iV 2 — 1) x 16 (iV 2 — 1) matrix defined by 

(^%„ = ^^ = -if.Ucr^. (20) 

Since the equation (|T9|) gives 10 (N 2 — 1) constraints among 10 (N 2 — 1) real variables, naively 
one would expect that only a very small portion of the configuration space survives. This is 
not the case, however. In fact, we find that all the configurations that belong to Qs satisfy 
the saddle-point equation (|I9D due to (p~8|) . The existence of solutions to flI9D other than of 
this type, cannot be excluded. It is reasonable, however, to consider that the configurations 
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in f2§ dominate in the sense of Lebesgue measure. Hence, at v = oo the full ten-dimensional 
Lorentz invariance is broken down at least to an eight- dimensional one. In order to examine 
whether the remaining eight-dimensional Lorentz invariance further breaks down, we still 
have to integrate over the saddle-point configurations. An important point here is that the 
imaginary part of the effective action for the saddle-point configurations takes only or 
7i due to (|13D . This means that there are actually two sequences of v, (I) v = 0, 2, 4, • • • , 2oo 
and (II) v — 1, 3, 5, • • • , (2oo + 1), which give two a priori different limiting theories. 

In order to formulate the integration over the saddle-point configurations, we first rotate 
them so that they satisfy Ag = Aiq = and then integrate over At, A2, ■ • ■ , A$. We define 
the Hesse matrix for r^[A] within the directions transverse to the integration domain as 

^=S=4^ ImlndetA1 =5 ImT^ ( A1 " e, ' , ° ,>, " e<M, ) ■ (21) 

where j, k = 9, 10. We find that the Hesse matrix Hj a ^b has 16 zero-eigenvalues with 
eigenvectors corresponding to perturbations 

(5A«, 5A a 10 ) = K, 0), (A a 2 , 0), • • • , (A a 8 , 0), (0, Af), (0, A 2 a ), • • • , (0, A a 8 ) . (22) 

These zero-modes are a reflection of the fact that the configurations after these perturbations 



still stay within f2 8 , and thus satisfy the saddle-point equation ( jT9|) . When N is even, the 
Hesse matrix Hj a ^ b has actually two more zero-eigenvalues. To see this, we first note that 
Hg a ,9b = Hi Qa:Wb = and H 9a:10b = H wb! Q a . Note also that H 9a:10b = —H 9btWa due to 
(|T0|). which means that an (iV 2 — 1) x (iV 2 — 1) matrix K^b defined by K ab = ifg^iofc is 
antisymmetric. Due to a general property of an antisymmetric matrix of odd size, K ab for 
even iV should have a zero-eigenvalue, whose corresponding eigenvector we denote as x a - 
Then one finds that (SAg, SA1 Q ) = (x a , 0), (0, x a ) are eigenvectors of the Hesse matrix Hj a kb 
with zero-eigenvalues. Unlike the 16 zero-modes in (P2"|), these two additional zero- modes, 
which exist only for even N, have nothing to do with the symmetry of the space of solutions 
to the saddle-point equation fllTf). Hence, they should be considered as accidental zero- 
modes. Since the saddle-point analysis in such a case becomes more complicated, we restrict 
ourselves to the odd iV case in what follows. 

In order to deal with the 16 zero-modes, which actually correspond to the Lorentz trans- 
formation in the (i, 9) and (i, 10) planes (i = 1, 2, ■ • ■ , 8), we have to take into account the 
phase volume analogous to the Faddeev-Popov determinant. In the present case, it can be 
given by the 16-dimensional volume spanned by the 16 vectors (^) in the configuration 
space. We denote this phase volume as V^ 16 ^. Thus we arrive at the models, which describe 
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the two limiting theories corresponding to the even/odd v sequences, 

V (X6)\j\-l/2 e -S b \ detM (S)\ ; (23) 



J v=2oo 




Z v=2oo+1 = / TTdAj |J|-^ e -^ detM^ , (24) 



where J denotes the determinant of the Hesse matrix Hj a ^b after removing the zero-modes. 
We recall that det is real but not positive definite, which makes the two models a priori 
different. One can see that the exact result ( |17D for (Sb) can be reproduced from (|23| ) and 
(H) as it should, by noting that V (16) i-> A 16 V (16) and J i-> \-2{2(v 2 -i)-i6} j under a gcale 
transformation A M i— > A/L^. 

Let us consider the question whether the remaining eight-dimensional Lorentz invariance 
of the models fl23|) and (24j) is further broken. For this we show that the hessian J becomes 
zero for configurations in and that the order of zero increases for configurations in 
with smaller d. To quantify this statement, let us consider a configuration with Ai, ■ • ■ , A^ 
being generic and Aa+i, ■ ■ ■ , A$ being of order e. We first note that each element of Hj a ^b for 
such a configuration becomes of order e^ 8- ^ due to (|T8|). Therefore, the non-zero eigenvalues 
of Hj a fib become of order e^~ d \ We have checked numerically that they are not of order 
higher than e 1 - 8- ^ generically. This means that J is of order e ( 8 - rf ) { 2 (^v 2 - 1 )- 16 }. On the 
other hand, degenerate configurations with smaller d is suppressed by the entropy factor 
e (s-d) (iv 2 -i) com j n g f rom the integration measure. The phase volume factor V*- 16 -* gives also 
a suppression of order e 2( - 8 ~ d \ Collecting all the powers of e, we find a suppression of 
order g 10 ( 8 ^ d ). This means, first of all, that the integrals (|23D and (24) are non-singular at 
the degenerate configurations. Secondly, we note that, as far as the iV-dependent part is 
concerned, the suppression for smaller d coming the entropy factor e 1 - 8 ^ ( Ar2_1 ) is exactly 
cancelled by the enhancement coming from \J\- 1/2 . This means that there is essentially no 
entropical barrier against having smaller d. 

So far, we have been focusing on the effect of on the dynamics of the IIB matrix 
model. The effect of Sb and r^, on the other hand, has been studied in Ref. || by using the 
low-energy effective theory obtained by the one-loop perturbation theory (||). There, it was 
argued that the complicated branched-polymer interaction among the diagonal elements Xj M 
of the bosonic matrices A^ might induce a collapse of the distribution of x^. Monte Carlo 
simulation of the low-energy effective theory at v = (TB| shows that the effect of <S& and 
is not sufficient to induce a spontaneous breakdown of Lorentz invariance. However, 
after taking into account the effect of by sending v to infinity, it is very plausible that 
the remaining eight-dimensional Lorentz invariance of the models ( p3|) and (p4|) is further 
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broken by the effect of S& and T^ r \ because there is no entropical barrier against having 
degenerate configurations any more. If this is the case, the vacuum at v = oo is given by 
degenerate configurations with d < 8. In this regard, we recall also that the pfaffian Pf M.(A) 
is zero when {A^} £ Q2, as stated in (15) . Therefore, the dimensionality d of the vacuum 
configurations must be d > 2. 

In order to check the above statements and to determine the dimensionality d of the 
vacuum configurations at v — 00, one has to carry out the integration over the saddle-point 
configurations described by (|23| ) and ([24]), for example, by Monte Carlo simulation. Note 
that the model (^) is not plagued by the sign problem any more. The model (|24]) , on the 
other hand, can be studied using the configurations generated with the model (|23|) as 

If the fermion determinant det does not have a definite sign for dominant configurations, 
the denominator as well as the numerator of the r.h.s. in fl2"5D becomes very small, exhibiting 
the sign problem. However, if it turns out that the vacuum of the model ( |2"5D is given by 
configurations with d < 6, then due to (|T4]), we have sgn(det .M®) = 1 for dominant 
configurations. In this case, the sign problem, which a priori exists in the model ([24]) , is 
solved in a sense dynamically, and the model (|24] ) is actually equivalent to the model (|23|). 
since {0) v=2oo+1 = (0) u=2rx> for any observables O. 



To summarize, we considered a deformation of the IIB matrix model by introducing an 
integer parameter v which couples to the imaginary part of the effective action induced 
by fermions. We studied the deformed model at v = 00, where the integration over the 
bosonic matrices is dominated by the configurations for which is stationary. First of 
all, there is still a huge configuration space left as the saddle-point configurations which 
we have to integrate over. Secondly, these saddle-point configurations have more than two 
shrunken directions. Thirdly, the more shrunken directions the configuration has, the more 
stationary rW becomes. This gives rise to an enhancement for configurations with more 
shrunken directions, and the effect was shown to cancel exactly the iV-dependent entropical 
barrier against having such configurations. An intriguing feature of this enhancement is that 
it occurs exactly when the configuration becomes a lower-dimensional hyperplane. This may 
be responsible for generating a flat space-time instead of a curved one or a fractal one, as a 
result of the dynamics of the IIB matrix model. We argued that the dimensionality d of the 
space-time generated dynamically in the deformed IIB matrix model atz/ = oois2<<i<8 
and most likely d < 8. We derived the models which describe the integration over the 
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saddle-point configurations for the two a priori different limiting theories corresponding to 
the even/odd v sequences. Remarkably, the model with the even v sequence is not plagued 
by the sign problem. One can therefore study the model by standard Monte Carlo simulation 
to extract the dimensionality d of the vacuum configurations. If this turns out to be d < 6, 
the model with the odd v sequence is actually equivalent to the model with the even v 
sequence. 

Whether the original model [y = 1) belongs to the same phase as v = oo is a nontrivial 
question, which is not accessible through standard Monte Carlo simulation due to the sign 
problem. We quote, however, an example from history, where an analogous approach was 
successful. It is the strong coupling limit in the lattice formulation of nonabelian gauge 
theories. Although one should send the bare coupling constant to zero in the continuum 
limit, confinement in these theories has been clearly demonstrated in the strong coupling 



limit [15| . The limit also provides a qualitative understanding of the spontaneous breakdown 
of chiral symmetry . The existence of such an approach is indeed one of the advantages of 
having a nonperturbative formulation, which should also apply to the case of string theory. 
Note, however, that the above successes in the case of nonabelian gauge theories rely crucially 
on the fact that there is no phase transition between the strong coupling regime and the 
weak coupling regime ||17|| . An illustrative example clarifying this point is that confinement 
holds in the strong coupling limit even for abelian gauge theories, for which there is actually 
a phase transition to a deconfining phase at an intermediate coupling constant. 

The fact that the model at v = oo still has a huge configuration space to integrate over, 
which is quite peculiar to this system, may suggest that it is in fact quite close to the original 
model {y — 1). Although it is hard to justify this statement rigorously, it would be certainly 
worthwhile to explore further the dynamics of the deformed IIB matrix model in this limit. 
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